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In contrast to finite dimensions where disordered systems display multifractal statistics only at 
criticality, the tree geometry induces multifractal statistics for disordered systems also off criticality. 
For the Anderson tight-binding localization model defined on a tree of branching ratio K = 2 
with A*' generations, we consider the Miller-Derrida scattering geometry [J. Stat. Phys. 75, 357 
(1994)], where an incoming wire is attached to the root of the tree, and where outcoming 
wires are attached to the leaves of the tree. In terms of the transmission amplitudes tj, the 
total Landauer transmission is T = \tj\'^ , so that each channel j is characterized by the weight 
Wj = \tjf/T. We numerically measure the typical multifractal singularity spectrum f{a) of these 
weights as a function of the disorder strength W and we obtain the following conclusions for its left- 
termination point a+{W). In the delocalized phase W < Wc, a+{W) is strictly positive a+{W) > 
and is associated with a moment index q+{W) > 1. At criticality, it vanishes a+{Wc) ~ and is 
associated with the moment index q+{Wc) = 1. In the localized phase W > Wc, a+(W) = is 
associated with some moment index q+{W) < 1. We discuss the similarities with the exact results 
concerning the multifractal properties of the Directed Polymer on the Cayley tree. 

I. INTRODUCTION 

Since its discovery fifty years ago Jj] Anderson localization has remained a very active field of research (see for 
instance the reviews 043)- According to the scaling theory Q, there is no delocalized phase in dimensions d = 1,2, 
whereas there exists a localization/delocalization at finite disorder in dimension d > 2. To get some insight into this 
type of transition, it is natural to consider Anderson localization on the Cayley tree which is expected to represent some 
mean-field limit. The tight-binding Anderson model on the Cayley tree has been thus studied by various techniques 
over the years [ol-fl^. Other studies have focused on random-scattering models on the Cayley tree [T5l-[T7|. For the 
version of the model defined on random regular graph of fixed degree, we refer to the recent work [l8| and references 
therein. The motivation to study Anderson localization on the Cayley tree has been revived recently by the question 
of many-body localization ^19"] , because the geometry of the Fock space of many-body states was argued to be similar 
to a Cayley tree ^ 26]. But of course, the questions on many-body localization are much more difficult and are still 
debated in the recent studies [22, [25l - l3l| . 

In quantum coherent problems, the most appropriate characterisation of transport properties consists in defining a 
scattering problem where the disordered sample is linked to incoming wires and outgoing wires and in studying the 
reflection and transmission coefficients. This scattering theory definition of transport, first introduced by Landauer 
[3^ . has been much used for one-dimensional systems [33l - (35| and has been generalized to higher dimensionalities 
and multi-probe measurements (see the review [3^ and references therein). For the Anderson model on the Cayley 
tree, an appropriate scattering geometry has been introduced by Miller and Derrida flS^ to perform weak-disorder 
expansions and numerical computations : an incoming wire is attached to the root, and outcoming wires are 
attached to the leaves of a tree of branching ratio K with N generations. In a previous work |14|, we have used this 
scattering geometry to study numerically the statistical properties of total the Landauer transmission T = \tj\'^ 
as a function of the number N of generations and of the disorder strength and to measure its critical behavior. The 
aim of this paper is to characterize the spatial inhomogeneity between the various channels j : the weights \tj\'^/T 
of the channels turn out to present a multifractal statistics, not only at criticality but also in the localized and 
delocalized phases as a consequence of the tree geometry. So we analyse how the singularity spectrum /(a) changes 
as a function of the disorder strength. 

The paper is organized as follows. In section [ill we introduce the Anderson localization tight-binding model on 
the Cayley tree and the scattering geometry that we consider to study the multifractal statistics of the Landauer 
transmission. Our numerical results concerning the multifractal statistics in various phases are described in sections 
mil and IIVI for the Box distribution and for the Cauchy distribution of disorder respectively. Our conclusions are 
summarized in section [V] In Appendix \^ we recall the exactly known results concerning the multifractality for 
the Directed Polymer on the Cayley tree, as a classical model which is useful to consider as a comparison, both for 
conceptual and numerical purposes. Appendix |B] explains how the numerical singularity spectra presented on figures 
have been obtained. 
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II. SCATTERING GEOMETRY FOR ANDERSON LOCALIZATION ON THE CAYLEY TREE 
A. Anderson tight-binding model on the Cayley tree 
We consider the Anderson tight-binding model 

i? = ^e.|i><i|+ J2 \i><j\ 

i <i,j> 



(1) 



where the hopping between nearest neighbors < i,j > is a constant V = 1 and where the on-site energies are 
independent random variables drawn from the 'Box' distribution 



I „f W W 

The parameter W thus represents the disorder strength. We have also studied the case of the Cauchy disorder 

W 

PCauchy 



7r(e2 + W^) 



(2) 



(3) 



B. Miller-Derrida scattering geometry 
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FIG. 1: Scattering geometry of Ref. [T^ : the disordered tree of branching ratio K = 2 starting at generation n = and ending 
at generation A'^ (on the Figure N = 4) is attached to one incoming wire and to outgoing wires. The total transmission 



is T^EJijf = 1 
outgoing wires. 



where r is the reflection amplitude of the incoming wire, and tj the transmission amplitudes of the 



We consider the scattering geometry introduced in [T3| and shown on Fig. [T] : the finite tree of branching ratio K 
is attached to one incoming wire at its root (generation n = 0) and to outgoing wires at generation N. One is 
interested into the eigenstate \^ > that satisfies the Schrodinger equation 



inside the disorder sample and in the wires where one requires the plane-wave forms 



(4) 



(5) 



These boundary conditions define the reflection amplitude r of the incoming wire and the transmission amplitudes tj 
of the j — 1, 2, ..K'^ outgoing wires. To satisfy the Schrodinger Equation of Eq. |3]within the wires with the forms of 
Eq. O one has the following relation between the energy E and the wave vector k 



E = 2 cos k 



(6) 
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To simplify the discussion, we will focus in this paper on the case of zero-energy E = and wave- vector k = n/2, 
because the zero-energy E = corresponds to the center of the band where the delocalization first appears when the 
strength W of the disorder is decreased from the strong disorder localized phase. From the conservation of energy, 
the total transmission T is related to the reflection coefficient |rp 

T^^lt.f = l-|rp (7) 

3 

We refer to fis'] for the results of a weak disorder expansion within this framework, and for a numerical Monte-Carlo 
approach to determine the mobility edge in the plane {E, W). In IJ] we have studied the statistical properties over the 
disordered samples of the total Landauer transmission Tjv at zero energy £^ = as a function of the disorder strength 
W and of the number N of generations. In the localized phase W > Wc, the typical transmission Tjf = e''^^" decays 
exponentially with the number N of generations 

N 



\n{Tr)^\nT^{W>W.) (8) 

where ^joc represents the localization length. In the delocalized phase, the typical transmission remains finite in the 
limit where the number of generations N diverges 



r!fP = e'"^~(^<^-^) ~ T^{W<Wc)>Q (9) 

The total Landauer transmission T is thus an appropriate order parameter of the localization transition at the mobility 
edge Wc- We refer to for more details on the critical behaviors of the localization length S^ioc and of the asymptotic 
value Too{W < Wc)- In the present paper, we wish to analyse the statistics of the contributions tj of the various 
channels to the total transmission of Eq. [7] as we now explain. 



C. Statistical properties of the weights of the outgoing channels 

In each disordered sample, we consider the weights 

and the 'analogs' of Inverse Participation Ratios 



(10) 



It is useful to introduce the multifractal formalism with respect to AI ~ (or equivalently the large deviation 
formalism with respect to the variable N = (In M)/{\nK)) : one defines the typical exponents T*'^P(q) as the exponents 
governing the decays of the typical values 



pyPiM = K^) = e''''^^''^ cx Af-'^^W = e-^(i"^)-''"'(9) (12) 

The typical singularity spectrum /*^^(a) is defined as follows : in a large disordered sample, the number Mm (a) of 
channels j (among the total of M of channels) that have a weight Wj scaling as Wj ~ M~" scales as 

A/-^P(a) ~ M-^""'^") (13) 
Then saddle-point computation of Ig (Eqs[Tl1and fT2)) 

M „ 

llyP{M = K^) = Y, - / dQM^"'"^")-^" (14) 
leads to the Legendre transform formula 

- T'yP{q) = max„ (a) - qa] (15) 
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Let us now briefly recall some basic notions about multifractality that will be useful to analyse the numerical 
results. As a consequence of the weight definition of Eq. [TUl the index a cannot be negative, so one has a> 0. As a 
consequence of Eq. [HI the 'typical' singularity spectrum is non- negative : py^{a) > 0. [Note that this is in contrast 
with the 'averaged' singularity spectrum /""(a) which can become negative to describe rare events (see fi\ for more 
details), but in this paper we only consider the typical singularity spectrum]. The terminations points a± are defined 
as the points where the singularity spectrum vanishes f{a±) — 0, whereas the singularity spectrum remains strictly 
positive in between 

fyP{a) > for a+<a<a_ (16) 

The left termination point a+ which represents the smallest possible a will play an essential role in the following. 
From the point of view of the Legendre transform formula of Eq. [151 it is associated with some positive value g+ > 0, 
where the saddle point a{q) reaches a+, so that for all higher q, the saddle point remains frozen at this value 

"(<? > 9+) = a+ (17) 

and the typical exponent t^^p^q) is simply 

T'yP{q >q+)^ qa+ (18) 

The same discussion can be transposed to the right termination point a_ associated with some negative index 
q„ < 0, with a{q < q^) — q;_ and T^yP{q < q ) ~ qa^. The value g = is associated with the most probable value 
ao = a{q — 0) where the singularity spectrum reaches its maximum 

/(ao) = 1 (19) 

Finally the value q = 1 is associated with the value ai = a{q ~ 1) where the singularity spectrum satisfies 

/(ai) = ai (20) 

as a consequence of the normalization I*^^ — 1 corresponding to T**'P(g = 1) = 0. 



D. Comparison with Anderson localization models in finite dimension 

We should stress here the similarities and differences with the usual multifractal definitions used for Anderson 
localization models in finite dimension d (see the review fl\) : for a normalized eigenfunction on the volume V — L''' 
(i.e. ^rev=L'' IV'('')P = l)i the Inverse Participation ratios are defined as 

p,^ I'^wi'' (21) 

reV=L'' 



and the exponents r(g) are defined as 



P'yP cx £--'""(9) (22) 



In finite dimension d, powers of L and powers of the volume V = L"^ correspond to the same scaling (up to a 
redefinition of the exponents), while on the tree, one should use powers of the number M — in the definitions of 
Eq. [T21 and not powers of the linear distance N . 

Another difficulty with the tree geometry is that sites of different generations are not equivalent in the pure case 
(see fl3| for explicit expressions of wavefunctions that decay exponentially with the distance N in the pure case). 
This shows that direct generalizations of Pq where the sum is over all sites of the tree is not appropriate (see again 
[T^ for more detailed discussion on the anomalous behavior of usual I.P.R.s), and this is why we have chosen to 
consider the weights of the channels in the Miller-Derrida geometry, since they involve the wave-function weights of 
the points that are at the same distance N of the origin. In the pure case, these weights have all the same weights 
|ijP/T = 1/ (see again [14] for more details), leading to the mono-fractal behavior of the Iq of Eq. [TT] 

|t|29 1 

T'i" = '<"') = /_S. .N. = lf«5riT (23) 



(ESiy^' 
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So the tree geometry has the pecuharity to induce multifractal behavior of the Iq even in the delocahzed phase (the 
radial symmetry of the pure case is not able to survive even at small disorder), whereas in any finite dimension, the 
I.P.R. in the delocalized phase are monofractal with the same scaling as the pure case. 

Finally in finite dimension, the localized phase is characterized by localized eigenfunctions where some rare sites have 
finite weights, whereas most sites have exponentially-small weights in the linear size L. On the tree, all eigenfunctions 
have to decay exponentially with the distance N , even in the pure case, to fullfill the normalization constraint with an 
exponentially-growing number of sites with the distance N . So in the localized phase, the fiuctuations of the weights 
will be also characterized by a multifractal statistics of the Iq of Eq. [TTl 

In summary, in contrast to finite dimensions where disordered systems display multifractal statistics only at criti- 
cality, the tree geometry induces multifractal statistics for disordered systems also outside criticality, if one consider 
the inhomogeneities among the points at a given distance from the center. In the recent mathematical study [37|. 
similar observables have been introduced with a large deviation analysis in N . 

Since these multifractal properties of disordered models defined on trees are unusual with respect to finite dimen- 
sions, it is useful to see how the multifractal analysis of Eq. [12] works in an exactly solved model : in the Appendix lAl 
we thus recall the case of the Directed Polymer on the Cayley tree, which is a classical disordered model having the 
same geometry. 



III. NUMERICAL RESULTS FOR THE BOX DISTRIBUTION 



In this section, we describe our numerical results for a tree of branching ratio K = 2 where the disordered on-site 
energies are drawn with the Box distribution of Eg. [21 The critical disorder width Wc at the center of the band E = 
has been found to be numerically in the interval (ol. [isj 

l6<Wc< 18 (24) 



A. Numerical details 



We have studied trees containing N generations with a corresponding number ns{N) of disordered samples with 
the values 

N = 10; 12; 14; 16; 18; 20; 22; 24 
n,{N) = 10^; 27.10^; 7.10^ 17.10^; 43.10^; 10*; 27.10^; 650 (25) 

We have chosen to work only with even N, because in the pure case, the total Landauer transmission is perfect 
rppure _ 2 QT^Yy for even N (see section 2.2 of Ref [l3| for more details). The transition amplitudes tj of the scattering 
eigenvalue problem of Eqs U and [S] are computed via the introduction of Riccati variables as explained in details in 
[isl [l^ . The multifractal spectrum is then obtained via the standard method of Ref 44 1, where the curve f{a) is 



obtained parametrically in the parameter q (see Appendix B [B] for more details) : here we have used values in the 
range —5 < q < +5. As shown in Appendix [K[ we have checked that the sizes and statistics of Eq. B5] were sufficient 
to obtain reliable results for the multifractal properties of the Directed Polymer model by a direct comparison with 
exactly known results in various phases. 

Let us make some final remark to explain the differences with respect to the numerical method used in our previous 
work concerning the full transmission T. In [l3| we had used the so called 'pool method' (see section 2.3.1 of Q|) 
which allows to study much bigger sizes (like N ~ 10^ generations). This was possible because the full transmission 
T can be directly computed from the refiection coefficient of the incoming wire alone (see Eq. (17) of [3); i-e. one 
does not need to compute the whole set of transmissions tj of individual channels. So the full transmission T can be 
obtained directly from the stable probability distribution of the complex Riccati variables, for which the pool method 
is well adapted. However the multifractal spectrum is a much more complicated observable : it does not depend 
only on the one-point distribution of the Riccati variables, but it involves the whole correlations between the Riccati 
variables along branches (each Riccati variable is computed from its K descendants, see Eq. (25) of [3])- Of course 
one could try to develop some generalized pool method to compute multifractal properties, but one should be very 
careful to avoid artifacts. In the present paper, we have thus chosen to work only with exact numerical data on finite 
trees to avoid any doubts on the numerical results. 
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FIG. 2: Box disorder W — 5 (a) The singularity spectrum f{a) has for termination points q+ ~ 0.37 and q_ ~ 2.82, for 
typical value ao — 1.33 and for tangent point ai — /(cn) — 0.76 (b) The corresponding a{q) saturates at a value around 
(7 I ~ 3. 
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FIG. 3: Box disorder W — 10 (a) The singularity spectrum f{a) has for termination points a+ ~ 0.1 and a- ~ 5.2, for typical 
value ao — 2.33 and for tangent point ai = /(cti) — 0.31 (b) The corresponding a{q) saturates at a value around q+ ~ 2. 



B. Delocalized phase 



In the delocalized phase W < Wc, we find that the left termination point is strictly positive Q!+(VF) > and 
is associated with a moment index q+{W) > 1. Two examples of our numerical data are shown on Fig. [2] and [3] 
corresponding to W = 5 and = 10 respectively. 



C. Critical point 



At criticality, the left termination point vanishes a+{Wc) = together with the tangent point ai = /(q?i) = 0, as 
shown on Fig. U corresponding to W = 17. The corresponding saddle-point a{q) saturates at the value (7-|_(VFc) ~ 1. 



D. Localized phase 



In the localized phase W > Wc, the vanishing left termination point Q!+(Wc) = is associated with some moment 
index q+{W) < 1, as shown on Fig. [5] and [5] corresponding to = 30 and W = 40 respectively. 
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a q 

FIG. 5: Box disorder W — 30 (a) The singularity spectrum f{a) has for termination points a+ = and q_ ~ 7.1, and for 
typical value ao ~ 4.0 (b) The corresponding a{q) saturates at the value g+ ~ 0.55. 




FIG. 6: Box disorder W — 40 (a) The singularity spectrum /(a) has for termination points a+ ~ and a_ ~ 7.35, and for 
typical value ao — 4.3. (b) The corresponding a{q) saturates at the value q+ ~ 0.45. 



8 



IV. NUMERICAL RESULTS FOR THE CAUCHY DISTRIBUTION 

In this section, we describe our numerical results for a tree of branching ratio K — 2 where the disordered on-site 
energies are drawn with the Cauchy distribution of Eq. |31 The numerical details are the same as in the section UlI Al 
The critical disorder width Wc at the center of the band E = has been previously found to be numerically in the 
interval @ 

3.6 <Wc< 4.4 (26) 

In many areas, the Cauchy distribution whose variance is infinite leads to anomalous results with respect to bounded 
distributions. In the context of Anderson localization, the Cauchy disorder is of course anomalous from the point 
of view of weak-disorder expansion which contains explicitly the variance of the disorder (see [35| and references 
therein). However from the point of view of Anderson localization transitions at finite disorder, we are not aware of 
any statement concerning its anomalous behaviors (except old conclusions concerning the absence of transition that 
have been shown to be false afterwards). On the theoretical side, the Cauchy distribution is well known to have 
many advantages : it is the only disorder distribution that leads to an exact and simple solution in one dimension 
(see [1^ and references therein), and that leads to an exact and simple solution for the density of states in any 
dimension [s^. For the Anderson localization on the Caylee tree that we consider in the present paper, it is also the 
only disorder distribution that leads to an exact and simple solution for the stationary distribution of the Riccati 
variables in the localized phase [l3| , the only remaining problem being that it is not known in the delocalized phase 
where the Riccati variables are complex [l^ . These theoretical advantages of the Cauchy distribution justify to study 
numerically its properties and to compare with the case of bounded distributions. In the following, we obtain that 
the results concerning the multifractal properties of the transmission for the Cauchy distribution are qualitatively the 
same as the results obtained in the previous section concerning the box distribution. 

A. Delocalized phase 

In the delocalized phase W < Wc, we find that the left termination point is strictly positive a+{W) > and 
is associated with a moment index q+(W) > 1. Two examples of our numerical data are shown on Fig. [7] and [8] 
corresponding to = 0.5 and W = 1 respectively. 

B. Critical point 

At criticality, the left termination point vanishes a+{Wc) = together with the tangent point ai = f{ai) = 0, as 
shown on Fig. [HI corresponding to = 4. The corresponding saddle-point a{q) saturates at the value q+{Wc) — 1. 

C. Localized phase 

In the localized phase W > Wc, the vanishing left termination point a+{Wc) = is associated with some moment 
index q+{W) < 1, as shown on Fig. [TOl and [TT] corresponding to = 6 and W = 10 respectively. 

V. CONCLUSIONS 

In this paper, we have studied the multifractal properties of the Landauer transmission for the Anderson localization 
tight-binding model on the Cayley tree within the Miller-Derrida scattering geometry. We have explained why, in 
contrast to finite dimensions where disordered systems display multifractal statistics only at criticality, the tree 
geometry induces multifractal statistics for disordered systems also off criticality. As an example, we have recalled in 
the Appendix the exact results concerning the Directed Polymer on the Cayley tree. We have presented numerical 
results for the typical multifractal singularity spectrum f{a) of the channels weights as a function of the disorder 
strength W, both the the Box distribution and the Cauchy distribution of disorder. Our main conclusion concerns 
the left-termination point a+{W). In the delocalized phase W < Wc, a+{W) is strictly positive a+{W) > and is 
associated with a moment index q+{W) > 1. At criticality, it vanishes a-i-(VFc) = and is associated with the moment 
index q+{Wc) = 1. In the localized phase W > Wc, a+{W) = is associated with some moment index q+{W) < 1. 
These properties of the delocalized and localized phases are thus qualitatively similar to the exact results concerning 
the Directed Polymer on the Cayley tree. 
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FIG. 7: Cauchy disorder W = 0.5 (a) The singularity spectrum /(a) has for termination points q+ ~ 0.37 and a_ ~ 5.3, for 
typical value qq — 1-5 and for tangent point ai = /(ai) — 0.74 (b) The corresponding a{q) saturates at a value around q+ ~ 3. 
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FIG. 8: Cauchy disorder W = 1 {a) The singularity spectrum f{a) has for termination points q+ ~ 0.21 and q_ ~ 6.63, for 
typical value ao — 2 and for tangent point ai — f{ai) ~ 0.55 (b) The corresponding a{q) saturates at a value around q+ ~ 2. 




typical value ao — 4.21 and for tangent point ai — f{ai) — (b) The corresponding a{q) saturates at the value q+ ~ 1. 
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FIG. 10: Cauchy disorder W — 6 (a.) The singularity spectrum /(a) has for termination points a+ — and q_ ~ 10.48, and 
for typical value ao — 4.68 (b) The corresponding a{q) saturates at the value g+ ~ 0.7. 




FIG. 11: Cauchy disorder VF = 10 (a) The singularity spectrum f{a) has for termination points a+ — and q_ ~ 11.06, and 
for typical value ao ~ 5.2 (b) The corresponding a{q) saturates at the value q+ ~ 0.5. 

Appendix A: Reminder on the Directed Polymer on the Cayley tree 

1. Reminder on the thermodynamics 

The Directed Polymer on a Cayley tree with disorder has been introduced in [sl] as a mean-field version of the 
Directed Polymer in a random medium . The model is defined by the partition function 

ZAT^^e-^^^'^) (Al) 
c 

where the configurations C are the paths of N steps on a Cayley tree with coordination number K. The energy 
E{C) of a path is the sum of the energies of the visited bonds. Each bond has a random energy drawn independently, 
for instance with the Gaussian distribution 



(A2) 
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This model presents many similarities [38|, |40| with the Random Energy Model, introduced by Derrida in the context 
of spin glasses [4l|. It presents a freezing transition at 

Tc - , ^ (A3) 

The free energy per step 4>{T) coincides with the annealed free energy above Tc and is completely frozen below [s^ 

I T 1 

2r ^ ^2lf " 2T 



4>{T) - 0a„n(r) = -Tlnif - — = - — for r > (A4) 



0(r) for T<T, (A5) 

2. Reminder on the finite weights statistics in the frozen phase 

The configurations weights in the partition function (Eq. lAll) arc defined as 



The moments 

A'" 

Y,{N)^Y.^l (A7) 

1=1 

have finite disorder-averages in the frozen phase /^(T) — T/Tc < 1 for values q > ^{T) = T/Tc fi^l 

"^-^^P^. with ,iT) = L 
The density giw) giving rise to these moments 



= r^.illn ' .:U with ^i{T) = ^ (A8) 



Yq = / dww'ig{w) (A9) 

reads [Hi 

and represents the averaged number of terms of weight w. This density is non-integrable as w — ?> 0, because in the 
limit N ^ oo, the number of terms of vanishing weights diverges. The normalization corresponds to 



Yq^i = / dwwg{w) = 1 (All) 







3. Reminder on multifractal properties of the weights 



In the non-frozen phase, the Yq of Eq. IA7I vanish with the number M — K of configurations as power-laws with 
non-trivial exponents for averaged and typical values, and it is convenient to introduce the multifractal formalism 
of Eqs [121 [HI HSl In the frozen phase, the finite asymptotic values obtained for q > fJ,{T) in Eq. IA9I correspond to 
Tq = 0, but it is nevertheless interesting to define the multifractal exponents of Eq. [T^lfor |g| < ii. The fact that the 
multifractal formalism is appropriate to describe the weights statistics for all values of T, and its exact computation 
is explained in detail in j45l - t52l |. Here we simply recall the main results, and we refer to |45|-|52| for more details and 
discussions. The main point is that the moments Yq can be rewritten as in terms of the partition functions at inverse 
temperatures /3 and \q\l3 

Yq{M = K^) = ^'-1 e-*?^^' ^ ^ ZN{\q\fi) ^ g-/3(|g|i=^„(|g|0)-9F„(/3)) (^p^i2) 
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so that the typical exponents of Eq. [12] are directly given in terms of the free-energy per step of Eq. IA5I 

r'y^{q)^^[W{W)~im] (A13) 

In the following, we thus quote the final results for the typical singularity spectrum f{a) in the various phases, and 
we present the numerical results obtained on trees of sizes 10 < iV < 22 to show that such sizes are sufhcient to 
obtain reliable results by comparison with the exact forms (more details on the numerical procedure can be found in 
Appendix [B]) . 

a. Non-Frozen Phase ^ = ^ > 1 

In the non-frozen phase, the left and right termination points read 

/ 1^^ 

a+ = 1 - - 

V 

^(l + i)' (A14) 
and the typical singularity spectrum is exactly Gaussian on the interval a+ < a < a_ where it exists 

/r>T.(") = - «+)("- - «) (A15) 

The terminating values a± are associated with the values q± = In the interval g_ = — /i < q < q+ = +H, the 
value a{q) dominating the saddle point calculation of Eq. [15] is simply linear in q 

a{q- =-^l<q<q+ = +^J)={l + ^- ^ (A16) 



The typical value corresponding to g = in Eq. IA16I 



ao = (^1 + ^ (A17) 

is the point where the singularity spectrum reaches its maximum f{ao) — 1 (Eq. I19p . Finally the value q — 1 where 
the singularity spectrum is tangent to the line a — f{a) (Eq. [20]) correspond to (Eq. IA16P 

«i = 1 - ^ = /(«i) (A18) 
The numerical results shown on Fig. [12] are in agreement with these expressions for ^ — 2. 

b. Critical point fi = -2!- = 1 

In the limit /i = ^ — > 1+, the above singularity spectrum has the following properties : the left terminating point 
a+ of Eq. IA14I vanishes 

(A19) 

together with the tangent point to the line a — f{a) (Eq. IA18[) 

ai = /(ai) = (A20) 
The corresponding numerical results for /i = 1 are shown on Fig. 1131 
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FIG. 12: Directed Polymer in the non-frozen phase at /i = T/Tc = 2 (see section lA 3 a|) : (a) Singularity spectrum /(a) 
: the terminating points are q+ = 0.25 and a- = 2.25, the typical value is ao — 1.25, the line a = f{a) is tangent at 
/(ai) = ai = 0.75 . (b) The saddle-point a{q) remains frozen at a+ for q > q+ — 2 and to a_ for q < q- — —2. 




FIG. 13: Directed Polymer at the critical point for jj. = T/Tc — 1 (see section lA 3 bh : (a) Singularity spectrum f{a) : the 
terminating points are a+ = and q_ = 4, the typical value is ao — 2, the line a — f{a) is tangent at the origin /(qi) = ai = 
. (b) The saddle-point a{q) remains frozen at q+ for q > q+ — 1 and to a- for q < q- — —1. 




FIG. 14: Directed Polymer in the frozen phase at /i = T/Tc — 0.5 (see section lA 3 ch : (a) Singularity spectrum f{a) : the 
terminating points are a+ = and a_ — 8, the typical value is ao = 4. (b) The saddle-point a{q) remains frozen at a+ for 
q > q+ = 0.5 and to a_ for q < q- — —0.5. 
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c. Frozen Phase /i = ^ < 1 

In the frozen phase, the left termination point is zero 

a+ = (A21) 

and the right termination point is 

4 

a_ = - (A22) 
M 

On the interval = < a < a- where it exists, the typical singularity spectrmri is again exactly Gaussian 

f^yPj,^{a) = ^a{a^ ~ a) (A23) 

The left terminating value a+ = is reached for q+ ~ ^ and the right terminating value a_ is reached for = — /i. 
In the interval q- = —fi <i q < q+ = the value a{q) dominating the saddle point calculation of Eq. [15] is again 
linear in q 

2 2 

a{q- = -/i < g < g+ = = 1 (A24) 

The typical value corresponding to g = in Eq. IA24l is 

ao - - (A25) 

Note that here the value g = 1 is already is the frozen region 1 > g+ = /i, so that the singularity spectrum is 
completely below the line a — f{a). The numerical results shown on Fig. [Tulare in agreement with these expressions 
for ^ = 0.5. 



Appendix B: Details concerning the numerical evaluation of f{a) 

In each sample with N generations, we compute the M ~ weights of Eq. IIOI from which we obtain immediately 
the parameters Iq of Eq. [TT] The corresponding typical exponents T^yP{q) of Eq. [Tulare then obtained by the following 
three parameters fit of the average over disordered samples 



In Iq{M) ~ -T^yP In M + cq ln(ln M) + cj, (Bl) 

i.e. T^yP is obtained as the coefficient of the leading linear term. The presence of the subleading term (cq) is important 
only in the regions where Tq nearly vanishes Tg ^ 0, whereas in the regions where Tq is not small, a direct linear fit 
could be acceptable and would give nearly the same numerical value for Tq. 

The typical multifractal spectrum f^^P^a) could in principle be obtained from T^yP{q) by some numerical procedure 
to perform the Legendre transform of Eq. [151 but this method has a lot of numerical drawbacks [44|. We have thus 
followed the standard method of Ref [l^l, with the simplification that we consider only boxes of size L = 1 in the 
notation of Ref [4^ . i.e. we do not use any coarse-graining with various box sizes, but we analyse instead the scaling 
with respect to the total number of boxes M — from our data obtained of various sizes N . The main idea of Ref 
(i^ l is to construct the following normalized g— measures from the initial measure defined by the weights of Eq. [TUj 

-5^^ - (B2) 



Of course q = 1 correspond to the initial measure w^^ = Wj. The denominator corresponds to Iq of Eq. [TT] It is 
then useful to introduce 



M 



Fq{M) = -Y,wf^ In 



M 



Aq{M) = -Y^w^;^hnw, (B3) 
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Fq{M) represents the Shannon entropy of the q— measure, whereas Aq{M) represents the averaged log of the initial 
weight Inwj with respect to the g— measure. From Eq. IB21 one obtains immediately the simple relation 

Fg{M) = qAg{M) + \nIg{M) (B4) 

which, after the division by the scaling factor (InM), exactly corresponds to the Legendre transform relation of Eq. 
[T51 Numerically, one only has to compute Iq and Aq , whereas Fq can be immediately obtained from them by Eq. IB3I 
The averages over the disordered samples of the two observables of Eq. IB3I can be analyzed by the following three 
parameters fits 

Fq{M) ~ fq In M + ci ln(ln M) + c\ 

Aq{M) ~ aq In M + C2 ln(ln Af ) + 4 (B5) 

The leading coefficients (a^, fq) then constitutes a parametric representation of the typical singularity spectrum f{a) 
as q varies. The presence of the subleading terms (ci,C2) are important only in the regions where fq nearly vanishes 
fq ~ 0, whereas in the regions where fq is not small, direct linear fits could be acceptable and would give nearly the 
same numerical values for (a,, fq). 

All the singularity spectra /(a) shown on the figures correspond to the parametric representation (a^, fq) obtained 
by the procedure just described. We have also presented our corresponding data for Uq to show the freezing phenomena 
in the parameter q. In Appendix A, we find that the singularity spectra obtained via this numerical analysis are in 
agreement with the available exact results. 
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